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ABSTRACT 

B1 ade-element/momentum performance prediction codes 
are routinely used for wind turbine design and 
analysis. A weakness of these codes is their 
inability to consistently predict peak power upon 
which the machine structural design and cost are 
strongly dependent. The purpose of this study was to 
compare post-stall airfoil characteristics synthesi- 
zation theory to a systematically acquired wind 
tunnel data set in which the effects of aspect ratio, 
airfoil thickness, and Reynolds number were investi- 
gated. The results of this comparison identified 
discrepancies between current theory and the wind 
tunnel data which could not be resolved. Other 
factors not previously Investigated may account for 
these discrepancies and have a significant effect on 
peak power prediction. 

1. INTRODUCTION 

The cost of a wind turbine is influenced by the 
designer's ability to predict peak power. Uncertain- 
ty in the performance prediction codes for predicting 
peak power forces the designer to be conservative and 
oversize the generator and other mechanical compo- 
nents. This practice leads to higher machine cost 
which in turn increases the resulting energy cost. 
The ability to predict peak power is dependent on an 
accurate description of the airfoil lift and drag 
characteristics immediately after stall. As a rotor 
approaches peak power, aerodynamic stall progressive- 
ly envelopes the blade from root to tip. At peak 
power the root end of the blade is deeply stalled and 
the tip region, where most of the power is being 
developed, is passing through the angle of attack 
range of 15 to 25 degrees. Peak power is most sensi- 
tive to the airfoil characteristics over this angle 
of attack range. 

Blade-element/momentum performance prediction codes 
provide the principle means for rotor design and 
analysis. With the use of reasonable airfoil data, 
these codes have been able to provide good estimates 
of rotor performance provided a substantial portion 
of the blade is not stalled. Past experience has 
shown that, for a fixed pitch rotor, when a large 
portion of the blade is stalled the use of two-dimen- 
sional (2-D) airfoil data under-predicts peak power 
and shows a rapid drop off thereafter relative to 
measured power. A detailed discussion of this 

problem is contained in References 1 and 2. The 
under-prediction of peak power has motivated the 
development of 3-D post stall airfoil data synthesi- 
zation routines (References 3 and 4). Of these two 
routines, that of Reference 3 is more widely used 
because of its user friendly nature. Without any 
user interaction the routine automatically provides a 
smooth curve fit from where the 2-D airfoil data 


leaves off up through the highest angle of attack 
experienced by the blade during stall. The approach 
of Reference 4 is more subjective and requires the 
user to provide discrete airfoil data in the light 
stall region. Both of these routines are based on 
past studies which show that after stall the airfoil 
characteristics of lift and drag are strongly 
governed by blade aspect ratio and are relatively 
independent of the airfoil geometry and Reynolds 
number. After stall, lower aspect ratios result in 
lower values of lift and drag coefficients due to 
greater flow leakage around the tip that relieves the 
pressure differential across the blade. 

Although the 3-D post-stall airfoil data synthesiza- 
tion routines provided the needed direction to help 
solve the peak and post-stall power prediction 
problem an adequate post-stall airfoil data base was 
not available at the time of their development. 
Consequently, these methods were tailored to predict 
peak power based on comparisons for a few machines. 
Applied to other machines these methods were suspect 
when predicted and measured peak power were not in 
agreement. For this reason, it was widely held that 
the 3-D post-stall synthesi zation routines needed 
further verification against a comprehensive post- 
stall airfoil data base. In an effort to satisfy 
this need, post-stall airfoil data were acquired in 
the Texas A&M University wind tunnel. 

2. POST STALL WIND TUNNEL TEST 

To establish the data base to better verify the post 
stall synthesi zation theory, nonrotating wing 
sections were tested in the Texas A&M 7 x 10 ft wind 
tunnel. Post stall performance characteristics were 
established as a function of aspect ratio, airfoil 
thickness, and Reynolds number. For this purpose the 
test models had a one foot chord and were of the NACA 
44XX family of airfoils. Blade aspect ratios of 6, 
9, and 12 and « were tested for airfoil thicknesses 
of 9, 12, 15, and 18 percent at Reynolds numbers of 
250,000, 500,000, 750,000 and 1,000,000. The results 
from this study showed that the post stall airfoil 
characteri sties of Cl and Cq are a strong 
function of aspect ratio as expected and a weak 
function of airfoil thickness. In the post stall 
region the influence of Reynolds number was 
negligible for the range tested. A report 
(Reference 5) containing the results from this study 
and can be obtained from Rocky Flats Wind Energy 
Research Center (WERC). 

Two plots extracted from this report which illustrate 
the influence of aspect ratio on Cl and Cq are 
shown in Figures 1 and 2. Figure 1 shows Cg versus 
angle of attack for the four aspect ratios. For 
discussion purposes the figure can be broken down 
into three distinct angle of attack regions. In the 
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Fig. 1. Effect of Aspect Ratio on Drag Coefficients 
of the NACA 4418 Airfoil at RN =0.25 x 10 « 


Fig. 2. Effect of Aspect Ratio on Lift Coefficients 
of the NACA 4418 Airfoil at RN =0.25 x 10 6 . 


TABLE 1. POST STALL C L , Cq EQUATIONS 


Equations of Reference 3 


C Qmax = 1,11 + 

0.018 AR 

0 (a =90' ) 

(i) 

Cq = Bi sin 2 a 

+ B 2 cosa 

(a=15°to90‘ ) 

(3) 

where: 

B » = %ax 
B 2 = Cq s - 

C tWx sin2 “s 




COSa s 



Cl = Ai sin2a 

+ A 2 C0S 2 a 
sina 

(a=15“to90'j 

(5) 

where: 

Ai - Bi/2 

*2 ■ (C Ls 

- C Dn,ax sin “s 

> sinae 

COSa s ) -S 

COS a s 



Equations Derived from Texas A&M Data 


Cq^ = 1.00 + 0.065 AR £ (a =90') (2) 

(0.9 + t/c } 


Cq = since + B 2 cosa (a=27. 5*to90° ) (4) 


where: 

Bl = C Dmax 

82 = C °s ' %ax sina s 

COSa s 


a 

CD 

16* 

0.100 

20* 

0.175 

25° 

0.275 

27.5° 

0.363 


Same as equation 5, however the use of the new 
C[u for evaluating Aj and A 2 results in a 
stronger AR influence. 


alternate: C[_ * C pcosa (a=30’to90*) (6) 

sina 


with: Cj_=0. 9 £ a=30 w 


AR - aspect ratio 

Cq - drag coefficient 

Ci - lift coefficient 

t/c - non-dimensional airfoil thickness 

a - angle of attack 

a s - angle cf attack at stall, usually 15* 
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first region up to 15 degrees, the wings are 
unstalled and for performance prediction purposes 
similar 2-D airfoil data for this region is normally 
acquired from one of many airfoil data catalogs for a 
desired Reynolds number. However, for this 3-D 
testing induced drag effects steadily increase with 
angle of attack with the magnitude being^ greater for 
lower aspect ratios. This subtle trend is masked by 
the accuracy of the data over this angle of attack 
range. The second angle of attack region from 15 to 
27.5 degrees is worthy of special attention since it 
has the dominant effect on peak power prediction. An 
important characteristic of this region is that 
aspect ratio effects cannot be discerned. In this 
region the higher induced drag associated with lower 
aspect ratio is neutralized by the lower pressure 
drag associated with lower aspect ratio. Aspect 
ratio effects are not seen to have a net influence on 
drag coefficient until the angle of attack exceeds 
27.5 degrees. Above this angle the pressure force 
across the wing increases as the tip relief 
diminishes with aspect ratio. 

For comparative purposes the post stall synthesiza- 
tion equations of Reference 3 were used to generate 
corresponding values of Cl and Cq. The calculat- 
ed values are compared to the test data over the 
angle of attack range of 20 to 90 degrees as shown by 
the large open symbols for aspect ratios of 6 and 
12. For both aspect ratios the calculated values 
fall well below measured values. In addition, the 

equations lack the strong sensitivity to aspect ratio 
as shown by the data above a 30 degree angle of 
attack. To better reflect the measured first order 
aspect ratio and second order airfoil thickness 
effects, equations 2 and 4 are presented for 
Cn and Cq, respectively. These equations 
wert x used to generate the shaded symbols and show 
good agreement with test results over the angle of 

attack range of 27.5 to 90 degrees. The new maximum 
drag coefficient equation has three times greater 
sensitivity to aspect ratio. In addition further 
aspect ratio sensitivity was acquired by using a 
first order (Sin) term in the drag coefficient 
equation. The new approximation (Equation 4) only 
applies over the angle of attack range of 27.5 to 
90 degrees versus 15 to 90 degrees used in the 
original approximation (Equation 3). The discrete 
values of drag shown in Table 1 for the angle of 

attack range of 16 to 27.5 degrees are for airfoils 
having thickness of 12 percent or greater. The wind 
tunnel test results show that for thinner airfoils a 
more rapid increase in drag can be anticipated over 
this local angle of attack range. A step jump in 
drag is needed to transition from 2-D data to 

discrete 3-D data. An increase to a value of 0.10 
appears reasonable in going from 15 to 16 degrees. 

A comparison of the calculated measured lift coeffi- 
cient versus angle of attack is shown in Figure 2 for 
the various aspect ratios. Prior to stall a reduc- 
tion in aspect ratio reduces the slope of the Cl 
curve. In blade-element/momentum type analyses this 
trend is approximated through the use of the tip loss 
factor. In the post stall region the measured lift 
coefficient undergoes a rapid initial drop off 
followed by a moderate reduction in Cl up until the 
angle of attack region of 60 degrees. An approxima- 
tion of this post stall behavior is obtained using 


Equation 5 from Reference 3. The equation provides a 
continuous drop off from peak C[_ to zero at 
90 degrees in a manner characteristic of the test 
data. The equation makes no attempt to approximate 
the initial drop and recovery the Cl experiences 
for the higher aspect ratio blades. However, this 
simplification does not appear to result in signifi- 
cant error. Calculated values of Cl are shown by 

the large open symbols using this equation with the 
values of Ai and A 2 based on the old C(i 

(Equation 1) as given in Reference 3. Calculated 

values of C L using the new Cn (Equation 2) 

to determine A* and A 2 are shown by the solid 

symbols. Stronger aspect ratio dependency and 
somewhat better agreement is achieved using the new 
approximation for Cq . This equation provides 
good results from peaK Cl throughout most of the 

post stall region which typically begins around 

15 degrees for most airfoils. 

Two interesting characteristics of these data in 
Figures 1 and 2 is that although the magnitude of 
Cl and Cq are a function of aspect ratio, the 

ratio of Cl/Cq for a given angle of attack is 

found to be independent of aspect ratio throughout 
the deep stall angle of attack range over 

30 degrees. In addition, Cl/Cq or more appro- 
priately L/D, passes exactly through a value of 1.0 
at an angle of attack of 45 degrees as would a flat 
plate. These characteristics are shown by the square 
symbols in Figure 3. For comparison flat plate 

theory is represented by the solid line as calculated 
by the equation L/D = cosa/sinot. The data points of 
L/D are seen to approximate the ratio of cosa/sina 

over the angle of attack range of 30 through 

60 degrees. Beyond this point the data points fall 
above flat plate theory due to the suction developed 
by the leading edge curvature. Based on this trend 
equation 5 can be replaced by the flat plate theory 
alternate of equation 6. Using equation 5, the ratio 
of L/D will be somewhat on the high side throughout 
the angle of attack range and will not quite satisfy 
the observation that L/D = 1 at an angle of attack of 
45 degrees as does equation 6. 



ANGLE OF ATTACK (DEG) 

Fig. 3 Effect of Aspect Ratio on Lift/Drag Ratio in 
Deep Stall 


37 



3. PEAK POWER PREDICTIONS 

Predicted performance was calculated for three fixed 
pitch wind turbines using Rocky Flats blade-element/ 
momentum analysis (PROPSH) which is described in 
Reference 2. Using both the post stall airfoil data 
syntheslzation method of Reference 3 and the modified 
equations based on the Texas A&M wind tunnel data 
predicted peak power was compared to measured peak 
power. The machines analyzed were the Jay Carter 25, 
ESI-54, and the MQD-0. In each case predicted peak 
power was calculated on the low side of measurements 
by 3 to 15 percent using the modified equations based 
on the Texas A&M wind tunnel data. Better agreement 
was achieved using the synthesization method of 
Reference 3. Calculated peak power relative to 
measured power varied from being 2 percent on the low 
side to 8 percent on the high side. Overall the 
modified equations based on the wind tunnel data 
predicted peak power 3 to 12 percent lower than the 
equations of Reference 3 which were developed by 
tailoring the equations to correlate with actual peak 
power measurements. 

The better agreement achieved using the equations of 
Reference 3 is attributed mainly to the higher L/D 
resulting in the angle of attack range of 15 to 
30 degrees. As previously indicated this portion of 
the post stall region has the greatest impact on peak 
power prediction. Another factor not to be overlook- 
ed is the accuracy of Ci associated with the 
2-D data tables used in the calculations. Any error 
in the Cl utilized translates Into a similar 
error in thecal ue of peak power. 

Speculation as to the difference in L/D resulting 
from the lift and drag equations of Reference 3 
versus the lower L/D resulting from the modified 
equations based on the Texas A&M wind tunnel data 
yields some potential explanations. The Texas A&M 
wind tunnel tests were nonrotational . 
Aerodynamicists over the years have speculated that 
for rotating wings centrifugal force effects are 
present in the boundary layer that result in spanwise 
flow toward the blade tip. Spanwise flow is thought 
to enhance performance by delaying separation of the 
boundary layer. Another source of spanwise flow is 
the spanwise suction gradient that results from the 
rotating blades local oncoming velocity being 
proportional to the blade radius. For both of the 
these causes, the spanwise flow is from blade root to 
blade tip. 

A second scenario that warrants serious consideration 
when trying to account for low peak power predictions 

deals with elastic twist of the blade. Prior to 

stall many airfoils have little or no nose down 

(negative) pitching moment. When stall occurs in the 
angle of attack range of 15 to 30 degrees, the moment 
coefficient rapidly diverges. The resulting large 
increase in nose down pitching moment can result in 
elastic twist toward feather which would enhance the 
power output and delay stall. The twist effect would 
become greater with windspeed as the stall enveloped 
the blade while propagated outward toward the tip. 
The occurrence of just a couple of degrees of elastic 
twist with increasing windspeed would substantially 
enhance peak power and post peak power. Elastic 
twist is a real consideration the manufacturer should 
be aware of because of its strong influence on 


measured post stall power. The degree of elastic 
twist can be expected to vary from one machine to 
another depending on the amount of blade torsional 
stiffness. 

4. CONCLUDING REMARKS 

The Texas A&M post stall wind tunnel test of the NACA 
44XX series airfoils as influenced by airfoil 
thickness, Reynolds number, and aspect ratio yielded 
the following findings. 

• Reynolds number was found to have an insignifi- 
cant influence on the airfoil characteristics 
over the range of 250,000 to 1,000,000. 

• Airfoil thickness was found to have little 
effect on post stall Cl although it does 
influence Ci which is normally acquired 
from a 2-D airfoil data catalog. Increased 
thickness did provide some reduction in drag 
throughout the post stall region. 

• Aspect ratio was the dominant influence on the 
post stall airfoil characteristics. However, 
this influence was not readily discerned over 
the angle of attack range of 15 to 30 degrees 
where the lower induced drag associated with 
high aspect ratios is being neutralized by a 
corresponding increase in pressure drag. Above 
30 degrees both drag coefficient and lift 
coefficient increased with aspect ratio such 
that the resulting ratio of L/D was independent 
of aspect ratio. 

Applying these test results the post stall synthesi- 
zation equations of Reference 3 were modified. Peak 
power predictions using the newly modified equations 
and those of Reference 3 showed that: 

• the modified equations underpredicted peak power 
by 3 to 15 percent 

• the equations of Reference 3, because of their 
more optimistic post stall L/D ratio, provided a 
better approximation of peak power. 

Potential causes of the differences between measured 
peak power and that predicted using the modified 
equations based on the Texas A&M data are: 

• the inability of nonrotating wind tunnel test to 
represent the influence of radial flow effects 

• elastic twists effects resulting from the 
divergence of the post stall moment coefficient. 
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